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ABSTRACT 

A prescribed salinity field is inserted into a 10-level primitive 
equation ocean circulation model. The model has been developed and is 
being improved in order to study large-scale thermal anomalies observed 
in the Central North Pacific Ocean by the North Pacific Experiment (NOR- 
PAX) , The salinity field, based on observations along 160 W longitude, 
is independent of longitude and time, and is smoothed in the north-south 
direction to remove small-scale variations, A new equation of state 
which is a function of temperature, salinity and depth is used to deter- 
mine density in the calculation of pressure from the hydrostatic equation. 
The parameterization of vertical mixing is changed to account for the 
stabilizing effect of salinity and the supercooling of the surface layers 
at high latitudes during winter. 

The addition of salinity induced changes in the currents at high 
latitudes when compared to the model without salinity. A Rossby wave 
was excited in the middle latitudes which produced transient changes. 

'The convective adjustment process enabled significantly colder water to 
overlie wanner water in the regions where salinity increases with depth. 
Supercooled surface temperatures , encountered during winter at high lati- 
tudes , are handled through convective mixing, but f^orther investigation 
into the dynamics of sub-grid scale vertical convection at near freezing 
temperatures is recommended. 
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I . INTRODUCTION 



The ocean and atmosphere behave as a coupled system v/ith anomalies 
in one medium inducing anomalies in the other. Bjerknes (1966, 1969, 

1972) and Namias (1959, 1969, 1972) have investigated this interaction by 
relating observed sea-surface temperature (SST) anomalies to anomalous 
behavior in the atmosphere. Other investigations have been conducted 
with numerical simulations to test the response of the atmosphere to a 
prescribed ocean (Rowntree, 1972), the response of the ocean to a pre- 
scribed atmosphere (Haney, 1974) or the simultaneous effects of a coupled 
system (Manabe , 1969; Bryan, 1969; Wetherald and Manabe , 1972). Through 
observations and numerical simulations a greater understanding of the 
air-sea system is sought with improved forecasting of the atmospheric and 
oceanic variability an ultimate goal. 

The numerical model used in this investigation is that of Haney (1974) 
with subsequent improvements involving seasonal atmospheric forcing, para- 
meterization of convective mixing and surface wind stirring (Haney and 
Davies, 1976) and nonlinear eddy viscosity (Haney and Wright, 1975). The 
goal of this investigation is to further improve this model by introducing 
salinity as a new variable. The model, based on the hydrostatic and 
Eoussinesq approximation, is a 10-level primitive equation model in a 
closed rectangular basin with uniform depth of 4 kilometers. The basin 
which extends from the equator to 65 ^N and spans 90^ of longitude is an 
idealized representation of the North Pacific Ocean. The horizontal re- 
solution of approximately 2® of latitude and 3® of longitude is provided 
by 33 uniformly spaced grid points in the north-south and east-west 
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directions. In the vertical, the ten levels are located at 10, 32.5, 

62.5, 102.5, 162.5, 462.5, 900, 1700 and 3100 meters below the surface. 
Wind and differential heating which drive the ocean circulation are de- 
rived from prescribed but time varying climatological values of solar 
radiation, cloudiness, surface air temperature, surface specific humid- 
ity and winds. Until the present study, the density in the model was 
assumed to be a linear function of temperature only. 

The model has been developed and is being improved in order to study 
large-scale thermal anomalies observed in the Central North Pacific 
Ocean by the North Pacific Experiment (NORPAX) . 

During a '* spin-up" phase, the model is integrated for approximately 
240 years of simulated time to produce a "model climatology". The result- 
ing circulation and temperature patterns adequately represent the large- 
scale north-south temperature variations and the major ocean gyres; 
however, the Alaskan gyre is not depicted. Once the model climatology 
has been established, temperature anomalies are injected into the model 
which is then integrated forward in time with or without anomalous atmo- 
spheric forcing. 'The evolution and migration of these anomalies can be 
related to their horizontal or vertical extent and atmospheric conditions 
to determine the controlling processes. Investigations by Hunt (1975) 
and Shiver (1977) detail the use of this model in such studies of large- 
scale ocean variability and the major conclusions of these works are pub- 
lished in Haney, Shiver and Hunt (1978). 

The motivation for this investigation is to detemine the effects 
of adding a prescribed salinity cross section to the model. The salinity 
cross section is steady state and independent of longitude. This new 
salinity field modifies the density field through a newly introduced 
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equation of state. The resulting modifications to the density field 
afreet the hydrostatic pressure, and hence the horizontal pressure gra- 
dient, as well as the vertical stability. 

This introduction of a prescribed salinity cross section is expected 
to improve the model in several ways. First, as noted above, the present 
model climatology (without salinity) does not depict an Alaskan gyre. 

Part of this model's deficiency may be due to the somewhat idealized 
form of the wind forcing at high latitudes; however, temperature varia- 
tions are small and salinity can play an important role in the density 
structure and hence, the geostrophic flow. Thus it is expected that 
salinity may improve the model's circulation at high latitudes. The de- 
sirability of adequately simulating the Alaskan gyre in this meteorologi- 
cally significant area is apparent. 

In his work using the model in which density depended on temperature 
only. Hunt C1975) noted that cold SST anomalies could not be maintained 
in the winter at high latitudes but were immediately convective ly mixed 
downward. Cold SST anomalies are a physically realizable condition dur- 
ing winter in high latitudes. Surface cooling occurs but the increasing 
salinity wirh depth provides stability and maintains the cold surface 
anomaly. However, free convection occurs and must be modeled whenever 
the surface layers become more dense than the underlying layers. The 
model previously determined stability solely from the temperature profile. 
If the temperatiara decreases with depth (positive profile) , the profile 
is stable and convection does not occur. On the other hand, an constable 
profile (produced for example by a cold surface anomaly) is instanta- 
neously adjusted to neutral thereby destroying any cold surface anomalies 
which produced the instability. In reality, density, and therefore 
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stability, is a complex function of temperature, salinity and depth and 
a negative temperature profile can be stabilized by a positive salinity 
profile. In this manner, cold SST anomalies can be formed and maintained 
and one of the goals of this study is to make the model more realistic 
in this regard. Similarly, but with less apparent significance, a warm 
anomaly can be maintained below the surface layers. The inclusion of 
salinity in the model can play an important role in the study of large- 
scale thermal anomalies by improving the climatological circulation pat- 
tern and altering the vertical stability, and hence, the convective 
adjustment parameterization. 
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II. DATA COMPILATION AND HANDLING 



The scilinity cross section used in the model is steady state and in- 
dependent of longitude. The steady state assumption was necessitated 
by the lack of data. A seasonal variation in surface salinity tied to 
the variation in evaporation is known to exist in several geographical 
areas. Other temporal variations in surface salinity exist at high lati- 
tudes where ice formation or melting alters the salinity profile. How- 
ever, the extent, both vertically and horizontally, to which these 
variations influence the open ocean are not known and can not be ade- 
quately modeled (Sverdrup / 1942). 

Allowing salinity to varv only latitudinally was dictated primarily 
by simplicity and economy; however, the dependence in Y alone is justi- 
fiable. The significant gradients in salinity are north-south and ap- 
preciable deviations from this pattern are found only at the boundaries 
where land rimoff and boundary currents become important (Figure 1) 

CRiley and Skirrow, 1975). Since the model is primarily concerned with 
open ocean circulation, the assimption of longitudinal independence in 
salinity seems reasonable. An anticipated result of the inclusion of 
salinity into the model is the formation of an improved and more realis- 
tic subarctic gyre. This gyre could be created by the geostrophic flow 
induced by a north-south salinity gradient. 

The salinity cross section prescribed in the model was derived from 
actual measurements from research vessels . The series , Oceanic Observa- 
tions of the Pacific: (Year) (Scripps Institution, 1960; 1963a and b; 
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1965a, b,c and d) compilsd by Scripps Institution of Oceanography was 
the source of the measurements . Salinity profiles were selected to ap- 
proximate a cross section along 160 longitude and at the model's grid 
point spacing of 2.0 3*^ latitude. Due to the sparseness of data in some 
areas, a common time of year could not be selected; however, 75% of the 
measurements lie between late spring and early fall. Because the Pacific 
Ocean narrows and shoals above 60 data in this area are limited in 
both horizontal and vertical extent. To overcome the data insufficien- 
cies the profile at 60 was duplicated for the grid points between 60 °N 
and the model's boundary which is located at 65°N. 

Figure 2 shows the latitude -height cross section, taken from the 
above data sources, plotted and analyzed at the model's grid poinrs . 

The salinity cross section (Figure 3) found in Intermediate Waters of the 
Pacific Ocean (Reid, 1965) was used to verify the representativeness of 
the compiled cross section. Major features such as the low salinity 
tongue penetrating equatorward from the high latitudes and the high 
salinity central water mass are adequately represented. Areas of sharp 
horizontal gradients at 35 and 45 °N are also reasonably depicted. 

Although the major feattires of the salinity cross section were well 
represented, a smoothing technique was employed to remove the small-scale 
irregularities and improve horizontal consistency . A repeated , three 
point weighted average given by 



S. , = .25 3. T + .5 S. , + .25 3 . . , 

j ,k 3-1/k j ,k 3+1/k 



( 1 ) 



was selected. Boundary values remained unchanged. This smoothing pro- 
cess was performed up to five times and the resulting cross sections 
analyzed. Four repetitions of this averaging scheme removed small-scale 
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irregularities (less than two grid lengths) but retained the major fea- 
tures of the salinity structure (Figure 4) . The central water mass re- 
mains intact and the low salinity tongue is still identifiable. Hori- 
zontal gradients have been smoothed but maintain their location and a 
portion of their intensity. Since the cross section was composed of 
actual soundings, vertical consistency was assured and vertical smoothing 
was not considered necessary. 
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III. EQUATION OF STATE AND CONVECTIVE ADJUSTMENT 



The density of sea water as a function of temperature, salinity and 
pressure is usually computed using the Knudsen formula. Although the 
Knudsen formula is quite accurate, it is generally considered too inef- 
ficient for inclusion in numerical models . This inefficiency arises due 
to the wide range of temperature and salinity for which the Knudsen for- 
mula holds , Since the actual variation of temperature and salinity in 
the open ocean is known to be limited, several researchers have developed 
empirical formulas which are accurate and efficient but apply only to 
open ocean temperature and salinity ranges. Two such empirical formulas 
were considered for inclusion in the model. 

Bryan and Cox (1972) solved an overdetermined, simultaneous system 
of equations which related departures from the mean of temperature and 
salinity to the departure from the mean density at a given depth. That 
is : 





- 2 , 2 „ 

OSi OS^ OT^ ... 




^Ik 




“ *^0k 




O's^ . . . 




^^2k 




"2 - pQk 


6t 

m 


6s 

m ^ 











m>n 



( 2 ) 



where: - (T,^ - 



T , p and S , are the mean values at a given depth. 
Ok Ok Ok 
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The density is given by; 



Xij^6t + + x^j^C6t)^ + x^j^(6s)^ + x^^5t6s = [p(T,S,Zj^)-p^j^]x 10^ 

(3) 

The constant coefficients, x^, were computed at each depth by Bryan and 
Cox using an iterative process . Acceptable accuracy was obtained with 
as few as three coefficients (n=3) and 50 temperature, salinity and den- 
sity combinations (m=50) . A significant improvement in accuracy was ob- 
tained with nine coefficients (n=9 , m=50) . Figure 5 illustrates the 
accuracy of this method as compared to the Knudsen formula. The errors 
are always largest near the surface where the widest range of tempera- 
tures and salinities exist. 

Another empirical formula for the equation of state of sea water was 
developed by Friedrich and Levitus (1972) . This formula is similar to 
the Bryan and Cox formula in that it fits a polynomial in temperature 
and salinity to computed values of density determined by the Knudsen for- 
mula; however, the coefficients are continuous, rather than discrete, 
functions of depth. 

A quadratic equation given by 

2 

aCT,S,Z) = ^ ^ 

is highly efficient but is not sufficiently accurate over the expected 
range of temperature and salinity. Further expansion of the polynomial 
to 

2 3 2 

aCS,T,Z) = + C^T + C^S + C^T + C^ST + C^T + C^ST (5) 
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yields acceptable accuracy while expansion to more terms does not im- 
prove accuracy appreciably . 

The coefficients are themselves quadratic polynomials in -Z (depth) 
which fit a smoothed curve to the pressure dependency of density. Due 
to the reduced range of temperature and salinity with depth, the coeffi- 
cients more accurately represent the pressure dependence at increased 
depth. Figure 6 shows the pressure dependence of the coefficients and 
the rms error of equation (5) with respect to the Knudsen formula. In- 
creased efficiency without a decrease in accuracy can be accomplished 
by rever-ting to equation (4) at depths of two kilometers and below. The 
values for the coefficients used in equations (4) and (5) are given in 
Tables I and II, respectively. Note that Z is in kilometers and defined 
positive downward. 

The Friedrich and Levitus equation of state was selected for inclu- 
sion in the model since computation of density at the model's irregular 
depths is facilitated by the continuous depth dependence. Another advan- 
tage over equation (3) of this formulation is that the equation of state 
would not have to be changed if the location of the levels in the model 
were ever changed. The explicit depth, temperature and salinity depen- 
dence of density as set forth by equation (5) aids in the remodeling of 
the convective adjustment process as discussed below. 

Without salinity the model had convectively adjusted whenever the 
profile was slightly positive, neutral or negative. The incorporation 
of salinity into the model enables a negative temperature profile to be 
convectively stable when coupled with a stabilizing salinity profile. To 
determine the relationship between depth, temperature, salinity and den- 
sity, the total derivative was expanded into its components. 
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( 6 ) 



dz “ 3z 3t 3z 3s 3z ' 



where Z increases upward. For stability, ^£.0/ temperature pro- 

file must satisfy 



i£ > f- i£ - i£i / i£ 

3z — 3z ~ 3s 3z^^^ 3 t 

(1) (2) (3) (4) (5) 



(7) 



■The right hand side of equation (7) gives the temperature lapse rate 
which renders the density profile neutral = 0) . Terms (1) and (4) 

are known from the temperature and salinity profiles computed in the 
model. Terms (2) , (3) and (5) are computed from the equation of state 
[equation (5)]. For efficiency, terms (2) and (3) can be simplified 
through scale analysis; however, all components of teicm (5) can be of 
similar magnitude and must therefore be retained. The expansion of the 
terms in equation (7) and the scale analysis is detailed in Appendix A. 
If the inequality is violated (as calculated using finite differencing 
between two levels) , then the temperature profile is convectively ad- 
justed. The adjustment is accomplished by satisfying two constraints: 

1) the heat present in the two layers remains constant; and 

2) the new temperature profile equals the neutral temperature 
lapse rate, the right hand side of equation (7) , plus a small 
additional stabilizing lapse rate. 

In this manner, cold water may overlie warm water provided the proper 
salinity str-ucture exists . 
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IV. RESULTS 



A '*spin up” of 246 years of simulated time without salinity (Haney 
et al . , 1978) was used as initial conditions. The salinity cross section 
and equation of state discussed above were inserted into the model and 
the integration continued for seven years. Streamlines and temperature 
contours at various depths for summer (mid-July) and winter (mid- January) 
of the model climatology, year 246 (without salinity) and the years 249 
through 253 (with salinity) were plotted, analyzed and compared. During 
the first seven year integration period, ’’neutral stability” was still 

3t 

defined as ~ 0- That is, the stabilizing effect of salinity in the 
parameterization of vertical convection (convective adjustment) was not 
included. This effect was added after year 253 (see below) . It was 
handled in this manner in order to study the effects of one change at a 
time. The model climatology, without salinity. Figures 7 through 11, 
depict the features previously discussed. In these and subsequent stream- 
line figures, the lengths of the small arrows are indicative of the rela- 
tive magnitude of the current at the grid point, while the long continuous 
lines are streamlines which help display the flow direction. In both the 
winter and sximmer the Equatorial Counter Current , North Equatorial Cur- 
rent, Kuroshio Current and North Pacific Current are well represented. 

At high latitudes a cyclonic turning is apparent but a closed Alaskan 
gyre is not formed. The temperature contours, Figures 12 through 16, 
illustrate the baroclinic zone formed by the North Pacific Current and 
the northward displacement of this baroclinic zone by the Kuroshio Current. 
The contour inteir/al is 3"^C at all levels except at level 9 where the 
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interval is 0,2^C. The North Pacific Central Water Mass is well defined 



in horizontal and vertical extent by the 21®C isotherm. 

After seven years of integration with salinity included, several 
changes in the streamline and temperature contours are visible. Compar- 
ing Figures 17 through 26 to their counterparts in the model climatology 
reveals a lessening of the cyclonic tendencies of the surface layers at 
high latitudes in the model with salinity included- Since this is 
directly opposed to the anticipated results , it merits further investi- 
gation. 

In the previous model with the equation of state represented by a 
linear function of temperature, the north-south gradients of temperature 
near O'^C, although small, generated significant geostrophic currents. 

Wind generated currents provided the cyclonic tendencies in the surface 
layers . The present equation of state is only weakly affected by the 
temperature fluctuations at low temperatures. The salinity gradients 
at high latitudes are very weak (recall, there is no horizontal salinity 
gradient from 60 to 65 "^N) with the resulting current primarily induced 
by the wind fields (Figure 27) . A strongly zonal forcing during winter 
generates meridional surface currents at high latitudes with the converse 
applying during simmer. Below the surface, zonal currents resulting 
from the north-south temperature and salinity gradients are encountered 
but their magnitudes are small. 

Another feature of interest, which appears during winter (Figures 
17a through 21a) when salinity is included, is the southward directed 
current along the western boundary near 45 ®N. This current is physically 
located in the region of the Oyashio Current, It appears only in winter, 
and since salinity is constant all year, atmospheric forcing is the 
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logical cause. The transport stream function (the vertically integrated 
mass transport is nondivergent and therefore represented by a stream 
function in this "rigid lid" model) for winter as ccmpared to summer 
(Figure 28) explains the seasonal dependence of this gyre. A strong cy- 
clonic stream function is located in this region in the winter while the 
value of the stream function in this area is small during summer. The 
salinity cross section (Figure 4) below 125 meters enhances the cyclonic 
shear due to a local salinity maximum coincident with the axis of the 
stream function miniinum. The resulting streamlines depict a geostrophic 
gyre below 100 meters and a wind induced gyre above 100 meters . When the 
cyclonic atmospheric forcing diminishes in the summer, the upper levels 
revert to a zonal flow and the gyre disappears except at the lower levels . 

A third change noted with the inclusion of salinity is a rather large 
cyclonic vortex in the eastern portion of the domain at about 35 ^N, A 
time series of streamlines from July 248 through July 251 (Figures 29 
through 32) show that the vortex is propagating westward at a rate of 
approximately 1.6 cm/sec. The vortex is present both winter and summer 
and penetrates from the surface to approximately 350 meters . ‘The charac- 
teristics of this vortex point to an internal Rossby wave. The phase 
speed of an internal Rossby wave is given by 



c 



u - 



2 . 2 

-.2 , 2 , o 

1 t k t :r- 
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( 8 ) 



where U = speed of the mean current 

1, k and m are the wave numbers in the x, y and z directions 
and 1 = 2TT/L 

f^ = average value for the Coriolis parameter 
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B = change in f with latitude 

3 t h 

N = Cag -^) the Brunt-Vaisala frequency 

This equation is based on a constant mean flow, U, and a Briint-Vaisala 
frequency, N, which is qualitatively applicable to the region in question. 
Values for the above equation were determined from the streamline plots , 
temperature profiles and velocity fields. The following values were 
selected: 

f = 2Q sincj) = 2f2 sin 35.5 = 3.4 x 10 ^ sec ^ 

0 

B = ~ ~ sec ^ cm ^ 

36 ^ 

lllxlO^cm 

1 = 2T/1000 km 
k = 2TT/3000 km 
m = 27r/300 m 

= 2x10"“^ X 980 X 8/2250 cm = 7x10”"^ sec~^ 

U = 1.5 cm/sec 

Inserting these values into equation (8) yields a value for the phase 
speed of -1.5 cm/sec which compares favorably with the obseirved speed. 

At this speed the wave would propagate across the domain of the model in 
12 years . 

The effects of the Rossby wave are well illustrated by the changes 
in the temperature contours (Figures 33 through 36). With the approach 
of the wave from the east, the baroclinic zone is displaced southward. 
While the baroclinic zone nearly returns to its original location after 
the wave passes, the temperature structure of the high latitudes is al- 
tered for a more extensive period. The high latitude isotherms are dis- 
placed southward and eastward; however the streamlines indicate flow from 
the south which should aid in returning the isotherms to their previous 
configuration . 
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The Rossby wave is therefore considered to be a transient feature 

due to the sudden insertion of salinity into the model and the salinity’s 

dynamical effects on the currents . As the wave slowly propagates across 

the basin, its amplitude decreases due to the lateral eddy diffusion in 

the model- The e-folding time scale for damping is of the order of 
2 2-1 

[ (k + 1 ) A ] , where A is the lateral eddy diffusion for heat. With 

ri H 

7 2 

A^ ^ 10 cm /sec, the time scale for damping is of the order of 5-10 years. 
Thus it appears that another 5-10 years of model integration will perhaps 
be necessary in order to allow the large-scale geostrophic currents to 
completely adjust to the new dynamical constraints imposed by the salinity 
field. This is not considered a serious problem. 

With the dynamical effects of the salinity reasonably well defined, 
the convective adjustment parameterization discussed in the prior section 
was inserted into the model. Within 40 days of simulated time the tempera- 
tures exceeded the F6.2 format. The individual terms of equation (7) were 
rechecked to insure their magnitudes and signs were correct. All terms 
produced appropriate values for realistic salinities, depths and tempera- 
tures; however, in the presence of strong surface cooling at high latitudes 

during the winter of year 253-254, unrealistic sea temperatures were 

9o 

generated. Temperatures below -3®C caused to become positive which 

3t 

forced the neutral value of to become large and positive. The para- 

oZ 

meterization of convection permitted temperatures below -3*^0; however, 
the equation of stats does not hold for temperatures below -3®C and does 
not model the change of phase from water to ice which occurs at subzero 
temperatures . 

While distilled water reaches maximum density before freezing, sea 
water with salinities greater than 25% exhibit increasing density with 
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# 



decreasing temperature throughout the cooling process; therefore, super-- 
cooled water continues to be convectively mixed (Sverdrup e^ al . , 1942) . 
Should the water column become supercooled through its entire depth, an 
ice cover would form insulating the column against further heat loss. 

The convective adjustment parameterization was altered to prevent 
the instabilities previously encountered. Water cooled below -2^C was 
convectively mixed so as to conserve heat but maintain minimum tempera- 
tures at or above -2°C. Should the entire column cool to -2"^C no further 
cooling would occur. With these changes inserted, the model was run 
through a simulated winter. Surface temperatures rapidly dropped to 
-2®C at the higher latitudes and were consistently cooler than the under- 
lying layers, the stratification being provided by the imposed salinity. 
Because low temperatures have little effect on density, changes in the 
circulation patterns were not noted. The inclusion of a convective ad- 
justment process dependent on depth and salinity as well as temperature 
allows cold anomalies to exist and more accurately models the physical 
(convective) processes obseir/ed in the high latitude ocean during winter. 
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V. CONCLUSIONS 



The inclusion of a prescribed salinity cross section into the model 
produced limited changes in the large-scale ocean circulation patterns . 
The most significant change was the appearance of a cyclonic gyre in the 
northwest portion of the domain. While its seasonal fluctuation at the 
surface is due to atmospheric forcing , the local salinity maximum below 
125 meters induced cyclonic shear at this latitude. No cyclonic gyre 
was formed in the region of the Gulf of Alaska presumably due to the 
reduced magnitude of the cyclonic atmospheric forcing in the northeast 
portion of the domain. The lack of boctom topography and of irregular 
boundaries may also impede the formation of an Alaskan gyre since all cy- 
clonic vorticity, which is put into that region by the wind, is free to 
propagate westward all the way to the western boundary. 

The propagation of an internal Rossby wave across the mid-latitudes 

altered circulation and isotherm patterns; however, these changes were 

of a nonpeimanent nature. The effect of dissipative forces were apparent 

after six years. Based on the prescribed value of the lateral eddy dif- 

7 2 

fusion coefficient of 5x10 cm /sec, the Rossby wave will be dissipated 
in approximately 10 years . 

The convective adjustment parameterization allows negative tempera- 
ture profiles to exist when balanced by a stabilizing salinity profile. 

To prevent supercooling of the surface layers, a simplistic mixing and 
ice formation parameterization was required. Although reasonable results 
were obtained, a more thorough investigation of the process of vertical 
convective mixing at freezing temperatures is recommended. 
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Several problems were encountered because the model exrends to 65 °N. 



Geographically , the Pacific Ocean above 60°N can not be realistically 
represented by a basin 4 kilometers deep with a width of 90® longitude. 
Salinity data above 60 ®N are sparse and highly variable since land runoff 
and ice formation or melting significantly alters the salinity profile 
there. For these reasons, 60 ®N would be a more appropriate northern 
boundary for the model. 

While the inclusion of salinity in the model did not generate all 
the desired results, it has improved the model in specific areas. The 
inclusion of salinity and the improved convective adjustment process 
are logical progressions in the foimiulation of realistic ocean models 
which can be effectively used in the study of large-scale ocean dynamics. 
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APPENDIX A 



DETAILS OF THE CONVECTIVE 
ADJUSTMENT PAPAMETERIZATION 



The empirical equation of state fomulated by Friedrich and Levitus 
(1972) is given by; 



a(S,T,Z) = c, + C,T + C S + C T^+ C^ST + + C,ST^ . 

1 2 3 4 5 6 7 



The values for the coefficients, which are polynomials in -Z (depth), 
are given in Table I. Note that Z is in kilometers and defined as posi- 
tive downward. Changing Z to positive upward and taking the partial 
derivative of sigma with respect to Z yields : 



3 ”- *‘'3, ^'"45 

3z 3z 3Z 3Z 'W 3z 



^"=6 3 ^^^7 2 



Inseirting the values for the coefficients and performing a scale analysis 
on the equation yields : 



i£ 

3z 



-5.1215 .IZ + .0363T - .0016ZT + .0086S 
2 

-.0006T - .0002ST + . . . 

-5.1215 (10 ^ gm/cm^ km) 



= -5.1215 



10 ^ (10 ^ gm/cm^ 



cm) . 



.0002SZ 



Likewise, the differentiation and scale analysis for sigma with respect 
to S is given by ; 

3a 2 

Is = C3 ^ C^T C,T 

= .8056 + .C086Z + .OOOIZ^ - .0030T - .0002ZT + . . . 

= .8056 (10~^ gro/cm^%); 
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however, the expression for the change in sigma with respect to T can 
not be simplified to a single value through scale analysis. Retaining 
all terms which may be significant over the expected range of tempera- 
ture, salinity and depth, the effect of temperature on density is given 
by : 

+ 2C^T + C.S + 3C.T^ + 2C_ST 
3T 2 4 o 6 7 

= .0498 + .0363Z + .0008Z^ - .0152T - .0013ZT - .0033 - .0002ZS 

2 

+ .OOOIT + .OOOIST. 

The stable lapse rate is given by: 

3t , 9S x y 3cr 

or 

^ > (5.1215x10 - .3056 ^) / {.K, + :< T + K 3 t .OOOIT (T-hS)) 

oZ — dZ 1 2 3 

where 

= .498 + .0363Z + .0008Z^ 

= -.0152 - .0013Z 
= -.003 - .00Q2Z 
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TABLE I. Coefficients for the Friedrich and Levitus equation 
of state, equation (5). C. = a. + b.Z + c.Z^ where 
Z is in kilometers and positive downward. 



i 


a. 

1 


b. 

1 


c . 
1 


1 


-7.2169x10"^ 


+5.1215x10“° 


-5.012x10“^ 


2 


+4.9762x10“^ 


-3.6349x10“^ 


+7.85 3x1 0“^ 


3 


+8.0560xl0”^ 


-8.5540xlo“^ 


+1.070xl0“g 


4 


-7.5911x10“^ 


+6.4295x10“^ 


-1.397xl0“g 


5 


-3.0063x10”^ 


+1.9 365xl0“g 


-3.899x10 g 


6 


+3.5137xl0~g 


-3.9740xl0“g 


-5.695xlo““ 


7 


+3.7297x10“^ 


-2 .8108x10“° 


+1.147xl0~ 



TABLE II . 


Coefficients for the Friedrich and Levitus equation 

of state, equation (4). C. = a, + b.Z + c.Z^ where 

i X 1 1 

Z is xn kilometers and posxtive downward- 


i 


a. 


b. 


c . 




X 


X 


X 




-2 


-0 


-2 


1 


-9.2163x10 2 


+5.1140x10 ^ 


-4.692x10 ^ 


2 


+4. 33l4xl0“^ 


-3.5685x10“^ 


+7.689x10“^ 


3 


+8.0640xl0“‘; 


-8.6826x10"’^ 


+1.433xl0“g 


4 


-6.2723x10“:: 


+5.1351x10“^ 


-1.246xl0“g 


5 


-2.7762x10“^ 


+1.7792x10 ^ 


-3.985xl0“° 
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Figure 1. Surface salinity contours of the Pacific Ocean (Sverdrup et al , ^ 1942). 
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Figure 2. Salinity cross section along 160°W with no horizontal 
smoothing . 



33 



Depth fm) 



Figure 3 



40 * 50 ^ 53^N 

3 3.8 32:2 



200 



400 - 



600 -^ 



800 - 



1000 

1000 



2000 -^ 



3000 - 



4000 



500 < 




Salinity cross section along 160 (Reid, 1965). 
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Salinity cross section along 160 °W smoothed four times. 
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Figure 5. Average percentage error in sigma for the Bryan and Cox 
equation of state where the Knudsen formula is taken as 
standard (Bryan and Cox, 1972). 




Figure 6. Coefficients of the Friedrich and Levitus equation of state 
as a function of depth and rms error of this equation with 
respect to the Knudsen formula (Friedrich and Levitus, 1972). 
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igure 7. Streamlines and velocity vectors for surface layer, year 246 
(without salinity) winter (top) and summer (bottom) . 
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Figure 8. Level 3 streamlines _ and velocity vectors 
year 246 (without salinity) winter (top) and 
summer (bottom) 
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Figure 9 . Level 5 
ear 246 (without 
umnier Cbottcm) 



streamlines and velocity vectors 
salinity) winter Ctop) and 
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Figure 10 . Level 7 streamlines and velocity vectors 
year 245 (without salinity) winter (top) and 
summer (bottom) 
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Figure 11. Level 9 streamlines and velocity vectors 
year 246 (without salinity) winter (top) and 
summer (bottom) 
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Figrore 12. Temperature contours for surface layer, year 246 Cwithout 
salinity 1 winter Ctop) and summer (bottom) . 
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Figure 13. Level 3 temperature contours, year 246 
(without salinity) winter Ctop) and summer (bottom) 








Figure 14. Level 5 temperature contours, year 246 
(without salinity) winter (top) and summer (bottom) 
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Figure 15. Level 7 temperature contours, year 246 
(without salinity) winter (top) and summer (bottom) 
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Figure 16. Level 9 temperature contours, year 246 
(without salinity) winter Ctop) and summer Cbottom) 








Figure 17. Surface streamlines and velocity vectors 
year 253 (with salinity) winter (top) and 
summer (bottom) 
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Figure 18. Level 3 streamlines and velocity vectors 
year 253 (with salinity) winter (top) and 
summer (bottom) 
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Fio-ure 19 . Level 5 streamlines and velocity vectors 
year 2 5.3 Cwith salinity) winter (top) and 
summer (bottom) 
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pio-ure 20. Level 7 streamlines and velocity vectors 
year 253 Cwith salinity) winter (top) and 
summer Cbottom) 
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Figure 21. Level 9 streamlines and velocity vectors 
year 253 Cwith salinity) winter Ctop) and 
Slimmer Cbottom) 
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Figure 22. Surface temperature contours, year 253 
(with salinity) winter (top) and summer (bottom) 
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Figure 23. Level 3 temperature contours, year 253 
(with salinity) winter (top) and summer Cbottom) 
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Figure 24 , Level 5 temperature contours, year 253 
(with salinity) winter (top) and summer (bottom) 
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Figure 25. Level 7 temperature contours, year 253 
Cwith salinity) winter (top) and summer (bottom) 
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Figure 26. Level 9 temperature contours, year 253 
(with salinity) winter (top) and summer (bottom) 
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Figure 27. Climatological surface wind stress 
patterns winter (top) and summer (bottom) 
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Figure 23. Transport streamfunction , winter Ctop) 
and summer Cbottom) 
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Figure 29. Level 3 streamlines and velocity vectors 
year 248 (with salinity) summer 
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30 . Level 3 streamlines and velocity vectors 
249 (with salinity) summer 




Q) 

tiO 0) 

■H 

CXH 



61 



31. Level 3 streamlines and velocity vectors 
250 (with salinity) summer 
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Figure 32. Level 3 streamlines and velocity vectors 
year 251 (with salinity) suimner 
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Figure 34. Level 3 temperature contours 
year 249 (with salinity) summer 
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Figure 35. Level 3 temperature contours 
year 250 (with salinity) summer 
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Figure 36. Level 3 temperature contours 
year 251 (with salinity) summer 







Figure 37. Climatological atmospheric quantities used in the 
model, (a) Solar radiation at the top of the atmosphere 
(Wm”'^); (b) fractional cloud cover (tenths); (c) surface air 

temperature (*C); Cd) surface vapor pressure (kPa) 
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